Did Swift measure GRB prompt emission radii? 
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ABSTRACT 

The Swift X-Ray Telescope often observes a rapidly decaying X-ray emission 
stretching to as long as t ~ 10 3 seconds after a conventional prompt phase. This 
component is most likely due to a prompt emission viewed at large observer 
angles 9 > 1/T, where 9 ~ 0.1 is a typical viewing angle of the jet and T > 100 
is the Lorentz factor of the flow during the prompt phase. This can be used to 
estimate the prompt emission radii, r em > 2tc/9 2 - 6 x 10 15 cm. These radii are 
much larger than is assumed within a framework of a fireball model. Such large 
emission radii can be reconciled with a fast variability, on time scales as short 
as milliseconds, if the emission is beamed in the bulk outflow frame, e.g. due to 
a random relativistic motion of "fundamental emitters". This may also offer a 
possible explanation for X-ray flares observed during early afterglows. 

Subject headings: gamma-rays: burster 

1. Introduction 

Recently launched Swift satellite (Gehrels et al. 2004) together with a network of ground 
based observations have been providing scientific community with crucial information on 
Gamma Ray Bursts (GRBs). Besides the landmark detection of afterglows from short GRBs 
(e.g. Gehrels et al. 2005), Swift has gathered crucial data on developments of GRBs at early 
times. This is especially important since early observations provide clues to the properties 
of the ejecta, like its composition, lateral distribution of energy etc. At late times the energy 
is mostly transfered to the forward shock, properties of which can hardly be used to probe 
the ejecta. A number of surprising results related to early afterglows have emerged (e.g. 
Tagliaferriet al. 2005; Nousek et al. 2005; Chincarini et al. 2005; O'Brien et al. 2005): (i) 
early, t < 10 3 s, rapidly-decaying X-ray component, (ii) X-ray flares occurring at t ~ 10 2 — 10 4 
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s, (iii) shallower than expected initial decay (or hump) of the afterglow. These features are 
common, but the light curves show a large variety. In this letter we discuss the first two 
mentioned effects, i.e. rapidly-decaying component and X-ray flares, since both can be related 
to the prompt emission (as opposed to afterglow) and can thus be used to probe the ejecta 
and the central engine. 

2. Prompt emission radii 

The initial fast-decaying part of afterglows can be a "high altitude" prompt emission, 
coming from angles 9 > 1/F (Kumar & Panaitescu 2000; Barthelmy et al. 2005), where 
9 is the angle between the line of sight and the direction from the center of the explosion 
towards an emitting point and V is the Lorentz factor of the outflow. For a 5-function in 
time prompt emission pulse, after an initial spike the observed flux should decay as t~( 2+Q ), 
where a ~ 0.5 is prompt emission's spectral index, (Fenimore et al. 1998), roughly consistent 
with observations. One also expects that prompt and early afterglow emission join smoothly, 
which seems to be generally observed (O'Brien et al. 2005). [Exceptions, like GRB050219a 
(Tagliaferriei al. 2005), may be due to interfering X-ray flares.] 

If we accept the interpretation of the fast decaying part as "high altitude" prompt 
emission, one can then determine radii of the prompt emission and compare them with 
model predictions. The currently most popular fireball model (e.g. Piran 2004) relates 
radii of emission r em to the variability time scale 5t of the central source r em ~ 2Y^c5t, 
where r ~ 100 — 300 is the initial Lorentz factor. Within the framework of the fireball 
model this is also the variability time scale of the prompt emission. Observationally, prompt 
emission shows variability on time scales as short as milliseconds, while most power is at a 
fraction of a second (Beloborodov et al. 1998). Adopting 5t ~ 0.1 s, the prompt emission 
radius is r em ~ 6 x 10 13 cm (r /100) 2 . If the emission is generated at r em and is coming to 
observer from large angles, 9 > 1/T, its delay with respect to the start of the prompt pulse 
is t ~ (r em /c)9 2 /2. If one can estimate 9, then this can be used to measure r em . This can 
be done from late "jet breaks", giving typically 9 ~ 0.1 (e.g. Frail et al. 2001). Then, for 
the X-ray tail of the prompt emission extending to t ~ 1000 seconds, the implied emission 
radius is r em > 6 x 10 15 cm. This is much larger than is assumed in the fireball model. To 
make it consistent with the fireball model and variability on short times scales, the Lorentz 
factor of the flow should be huge, r > 1000, but this would imply that emission is strongly 
de-boosted, T 9 ~ 100. Increasing 9 cannot save the day either since the required viewing 
angle would be 9 ~ 1, implying a jet moving always from an observer. 

Along the similar lines of reasoning, Lazzati & Begelman (2005) estimated prompt 
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emission radii for a particular case of GRB 050315 for which a possible jet break is identified 
(Vaughan et al. 2005) . Steep decay in that case is relatively short and lasts for 100 s, giving 
r em > 2.5 x 10 14 cm. Note, that any observed duration of the steep decay phase provides 
only a lower limit on the prompt emission radius since the end of the steep decay may be 
related to emergent afterglow emission and not to the fact that the edge of the jet becomes 
visible, see Fig. 1. On the other hand, late jet break time provides an estimate of the total 
opening angle of the jet. In any case, GRBs with longer lasting steep decay phase, up to 10 3 
s, provide the most severe constraints on the models. 

Thus, the interpretation of the fast-decaying initial X-ray light curve as prompt emission 
seen at large angles can hardly be inconsistent with the fireball model. We should then 
either look for alternative possibilities to produce the fast-decaying part of the X-ray light 
curve (e.g. Meszaros & Rees 2001), or consider models that advocate production of prompt 
emission at larger radii, see §5. 



3. Fast variability from large radii 

If prompt emission is produced at distances ~ 10 15 — 10 16 cm, how can a fast variability, 
on times scales as short as milliseconds, be achieved? One possibility, is that emission 
is beamed in the outflow frame, for example due to a relativistic motion of (using pulsar 
physics parlance) "fundamental emitters" (Lyutikov & Blandford 2003). To prove this point, 
we consider an spherical outflow expanding with a bulk Lorentz factor T with iV randomly 
distributed emitters moving with respect to the shell rest frame with a typical Lorentz factor 
7t- Highly boosted emitters, moving towards an observer, have a Lorentz factor 7 ~ 277T 
in the observer frame. If emission is generated at distances r em , the observed variability 
time scale can be as short as ~ (r em / c) /2^ 2 « (r em /c)/8(7rr) 2 , so that modest values of 
7t ~ 5 — 10 <C T ~ 100 — 300 would suffice to produce a short time scale variability from 
large distances r em ~ 10 15 — 10 16 cm. 

The model should satisfy a number of constraints. First, the number of sub-jets directed 
towards an observer from viewing angles 9 < 1/T should be larger than unity (in order to 
produce at least one true prompt emission spike), but should not be too large, otherwise 
prompt emission will be a smooth envelope of overlapping spikes. If a typical jet opening 
angle is 6j, then the number of sub-jets seen "head-on" from angles < 1/T is 

n promp t ~ (r7T 0. )2 - W 



This should be larger than 1. 
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The second constraint that the model should satisfy relates to the efficiency of energy 
conversion. Suppose that the thickness of an outflowing shell in its rest frame is L^eii ~ t s cT, 
where t s is a source activity time (t s ~ 30 — 100 s for long bursts and t s ~ 1 s for short 
bursts). Suppose then that fundamental emitters operate for a time t pu i se = f]tL s h e u in the 
flow frame, where r) t is a dimensionless parameter. During this time the source can tap into 
energy contained within volume (ct pu i se ) 3 . The ratio of this volume times the number of 
emitters to the total volume of the shell is a measure of efficiency of energy conversion into 
radiation: 

_ N(ct pu i 

V = r 2 OH cT [ ' 

1 em u j L s^ L 

Since tapping of energy in the volume (ct pu i se ) 3 is a definite upper limit on conversion effi- 
ciency, in the calculations we allow rj defined above to be slightly larger than unity. 

To produce light curves we calculate the intensity of emission from sub-jets that are 
randomly located within the shell and moving in random direction with random Lorentz 
factors 1 < 7r < 7r,max — 5. Each emitter is isotropic in its rest frame and is active for 
a random time < t' em < f]Tt s cT = t vu \se,max with t]t = 0.5. The observed intensity of 
emission from each sub-jet oc 5 3+a (Lind & Blandford 1985), where 5 = 1/7(1 — (3cos6 S j) 
is a total Doppler factor including bulk and random motion, 9 S j is an angle between the 
line of sight and direction of the sub-jet motion. As the burst progresses, larger angles and 
more of sub-jets producing prompt emission become visible. Most of them will be seen from 
angles > 1/tt in the bulk frame, producing a combined smooth curve overlaid with spikes. 
The average Doppler factor decreases with time 5 t s T/t and the average flux decays as 
£-(2+a) £-2.5 f or a _ 5 j n pjg^ 1 we pi Q t an example of a prompt light curve in this 

model. 



3.1. Lateral dependence of prompt emission 

Variations of the decay rate from the t"( 2+a ) law may be used to probe angular de- 
pendence L(9 ax i S ) of the intensity of the prompt emission, where d ax is is an angle between 
the axis of the explosion and an emitting point. More shallow decays can be due to, e.g. 
, a structured jet, with L ~ 0~^ is observed outside of some core: late time emission then 
is coming from the more energetic core part. The effective emission intensity increases ap- 
proximately as 9 2 oc t, and will result in an observed decay t~( 1+a K Similarly, if the prompt 
emission is seen within a core, late emission comes from less energetic wings, giving in case 
of a structured jet a flux oc t~( 3+a K Qualitatively, the relativistic internal motion of emitters 
makes it "easier" to see the high altitude emission. 



- 5 - 



To show this numerically we parameterize the number density of emitters as n(9) oc 
l/(9 2 + Oq), where 9 is an angular core radius. [There are, naturally other possible param- 
eterizations, e.g. of intensity of each emitter]. The results are presented in Fig. 2. 

We can also expect deviations from a simple power-law decay due to not exactly spherical 
form of the emitting surface. Such distortions are expected due to a development of the 
Kelvin-Helmholtz instability during an accelerating phase of the outflow. They won't be 
erased during the coasting stage due to causal disconnection of the flow separated by angles 
> 1/r. Additional complications may come from the way the data analysis is performed, 
e.g. through a choice of initial time trigger (Zhang et al. (2005), see also Lazzati & Begelman 
(2005)). 



4. Origin of X-ray flares 

Early X-ray light curves show complex behavior with flares and frequent changes in a 
temporal slope (e.g. O'Brien et al. 2005). Flares show very short rise and fall times, much 
shorter than observation time after the on-set of a GRB, while the underlying afterglow 
has the same behavior before and after the flare (Burrows et al. 2005) (though there are 
exceptions). Both of these observations argue against a physical process in the forward 
shock. In addition, there is a hardening of the spectrum during X-ray flares (Burrows et al. 
2005). 

In the present model we interpret X-ray flares as been due to sub-jets located at large 
viewing angles, 9 > 1/r, but directed towards an observer. Randomly located, narrow 
spikes are clearly seen in the model light curves, Figs. 1-2. In addition, as the flares are 
less de-boosted than the average high altitude outflow, they will have a harder spectrum, as 
observed. 



5. Discussion 

In this letter we first point out that the interpretation of the initial fast-decaying part of 
the X-ray GRB light curves as a prompt emission seen at large angles, and a generic estimate 
of jets' opening angle allows a measurement of the radius of prompt emission, which turns 

out to be relatively large, > 10 15 cm. On basic grounds, 7-ray emission should be generated 

/ \ 1/3 

before the deceleration radius r dec ~ ( ^pc 2 ™' 2 ) ~ ~~ cm ' wnen m ost energy of 
the outflow is given to the surrounding medium (here E iso is isotropic equivalent energy, p 
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is density of external medium, Y dec is Lorentz factor at r dec ). 1 The inferred emission radius 
is within this limit. 

The estimate of the emission radius is very simple, and, in some sense, generic. It can 
hardly be consistent with the fireball model, unless extreme assumptions are made about 
the parameters {e.g. very large Lorentz factor). On the other hand, there are alternative 
models {e.g. the electromagnetic model (Lyutikov 2005b), see also Thompson (2005)) that 
place prompt emission radii at large distances, just before the deceleration radius r dec . 

Secondly, we show how models placing emission at large radii may be able to reproduce 
a short time scale variability of the prompt emission and explain later X-ray flares. This can 
be achieved if the prompt emission is beamed in the rest frame of the outflow, which may 
be due to an internal relativistic motion of "fundamental emitters". 

What can produce a relativistic motion in the bulk frame? It can be due, for example, 
to a relativistic Burgers-type turbulence (a collection of randomly directed shock waves). It 
is not clear how such turbulence may be generated. Alternatively, relativistic internal sub- 
jets can result from reconnection occurring in highly magnetized plasma with a 1, where 
o is a plasma magnetization parameter (Kennel & Coroniti 1984). In this case the matter 
outflowing from a reconnection layer reaches relativistic speeds with 7 out ~ a (Lyutikov 
& Uzdensky 2003). Internal synchrotron emission by such jets, or Compton scattering of 
ambient photons, will be strongly beamed in the frame of the outflow. Note, that this model 
does not require late engine activity to produce flares. 

One of the main observational complications is that at observer times larger than the 
conventional prompt phase, the X-ray light curve is a sum of the tail of the prompt emission, 
coming presumably from internal dissipation in the ejecta, and the forward shock emission. 
It is not obvious how to separate the two components. For example, GRBs which do not 
show a fast initial decay may be dominated by the forward shock emission from early on 
(O'Brien et al. 2005). This uncertainty also affects estimates of the emission radius since 
the end of the steep decay may be related to the emergent afterglow emission and not to 
the jet opening angle (or observer's angle, in case of a structured jet), see Fig. 1. Another 
complication is that at these intermediate times, 10 3 < t < 10 4 s, even the forward shock 
emission itself often does not conform to the standard afterglow models, showing flatter than 
expected profiles {e.g. Nousek et al. 2005). 

A consequence of the model is that some short GRBs may be just a single spike directed 



1 Note, that rdec defined above is independent of ejecta content, contrary to the claim in Zhang & Kobayashi 
(2005), see Lyutikov (2005a). 
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towards an observer of a long GRBs. In our model the shorter spikes are highly beamed, less 
frequent and produce harder emission. This can apply only to some short GRBs since as a 
class they are well established to have different origin than long GRBs (from non-observation 
of a supernova signature and coming from a distinctly different host galaxy population). 

REFERENCES 

Barthelmy, S. D., et a/.2005, astro-ph/0511576 

Beloborodov, A. M., Stern, B. E., Svensson, R., 1998, ApJ, 508, 25 

Burrows, D. N. et al, 2005, astro-ph/0511039 

Chincarini, G. et al, 2005, astro-ph/0511107 

Gehrels, N. et al, 2004, ApJ, 611, 1005 

Fenimore, E. E., Madras, C. D., Nayakshin, S., 1998, ApJ, 473, 998 
Frail, D. A. et al, 2001, ApJ, 562, 55 

Gehrels, N., Barbier, L., Barthelmy, S.D., et al. 2005, Nature 437, 851 

Kennel, C. F., & Coroniti, F. V. 1984a, ApJ, 283, 694 

Kumar, P., & Panaitescu, A. 2000, ApJ Lett., 541, 51 

Lazzati, D., Begelman, M., 2005, astro-ph/0511658 

Lind, K. R., Blandford, R. D., 1985, ApJ, 295, 358 

Lyutikov M., Blandford R., 2003, astro-ph/0312347 

Lyutikov, M., Uzdensky D., 2003, ApJ, 589, 893 

Lyutikov, M., 2005, astro-ph/0503505 

Lyutikov, M., 2005, astro-ph/05 12342 

Meszaros, P. & Rees, M.J., 2001, ApJ Lett., 556, 37 

Nousek, J. A. et al, 2005, astro-ph/0508332 

O'Brien, P. T. et al, 2005, astro-ph/0601125 



- 8- 



Piran, T., 2004, Reviews of Modern Physics, 76, 1143 
Tagliaferri, G. et ai, 2005, Nature, 436, 985 
Thompson, C, 2005, astro-ph/0507387 
Vaughan, S. et al, 2005, astro-ph/0510677 
Zhang, B., Kobayashi, A., 2005, ApJ, 628, 315 
Zhang, B. et ai, 2005, astro-ph/0508321 



This preprint was prepared with the AAS IATgX macros v5.2. 



- 9- 



10 1 



10 10 et 



10 8 



10 8 r 



10 7 r 



10 6 



10 5 - 



10 4 



1000 



~~n r 



GRB phase 




X-ray flares 



? afterglow emission , 



R em /c)0 2 /2 



_i i i i i i i i i 1 1 1 1 i i i i i 1 1 1 




J I L-t-LJJ 



10 



100 



1000 

seconds 



10 4 



10 5 



Fig. 1. — Prompt emission produced by emitters moving randomly in the bulk frame. Emis- 
sion is generated within a shell of thickness t s c = 3 x 10 12 cm in observer frame, moving 
with T = 100 at distance r em = T 2 t s c by randomly distributed sub-jets with random orien- 
tation moving with random Lorentz factors 1 < 77- < iT.max — 5. Each emitter is active for 
random time < t' em < 0.5t s cr = t pu i seymax in its rest frame. Dotted line: average intensity 
oc t~( 2+a ) = t~ 2 - 5 . Dashed line: expected afterglow signal rising oc t 2 , peaking at ~ 100 s and 
falling off oc t~ L5 with arbitrary normalization. Homogeneous jet centered on an observer 
with opening angle 9 = 0.1, dimensionless parameters are n prompt = 1.2 and r\ = 1.6. 
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Fig. 2. — Same as Fig. 1 for structured jet with number density of sub-jets oc l/(9 2 + 9%), 
core size 9 = 7r/100, observer angle 9 b = tt/10, dimensionless parameters are n promp t = 1.26, 
rj = l.Q (calculated with 9 = 9 ob ). 



